CURRENT SHEET MAGNETIC MODEL FOR THE SOLAR CORONA Kenneth H. Schatten 


ABSTRACT A new magnetic model is developed and compared with previous models and the observed 
solar corona. An attempt is made to more accurately compute the three-dimensional 
currents flowing in the solar corona. Physical reasons are given that require most of the 
large-scale currents flowing in the solar corona to lie near thin sheets. The current sheets 
are not constrained into any particular geometry or symmetry as in the previous models 
of Altschuler and Newkirk [1969] and Schatten et al [1969] . A comparison with the 
axisymmetric, isothermal MHD solution of Pneuman and Kopp [1970] suggests that the 
model is able to simulate to high accuracy an isothermal corona. A comparison of the 
model with the May 30, 1965, solar eclipse and the November 12, 1966, solar eclipse 
shows the model is capable of computing many features including the polar plume orien- 
tations as well as radial and nonradial streamers in the solar corona. 


INTRODUCTION 

The advent of large digital computers and detailed mag- 
netograms has permitted sophisticated analyses of 
magnetic field configurations in the vicinity of the sun as 
suggested by Gold [1956] . Computations of the coronal 
magnetic field utilizing potential theory began with the 
Schmidt [1964] program to plot current-free magnetic 
fields above active regions. Rust [1966] has compared 
the field configuration of the Schmidt program with 
direct observations of prominent material. Newkirk et al 
[1968] utilized potential theory over the entire sun to 
calculate field patterns for a comparison with the pro- 
jected appearance of the November 12, 1966 solar 
eclipse. Schatten [1968 a,b] and Schatten et al [1969] 
developed a “source surface” technique to calculate the 
effect of coronal currents upon the field. The currents 
were chosen to draw the field into a radial direction 
[fig. 1). This model allowed comparisons of fields cal- 
culated with the interplanetary field, a Faraday rotation 
eclipse [, Schatten , 1968 a,b 1969, 1970; Stelzried et al, 
1970; and Smith and Schatten, 1970] . The technique has 
received favorable review by Cowling [1969] . 

Schatten et al [1969] utilized a “source surface” 


The author is at the Laboratory for Extraterrestrial Physics 
NASA-Goddard Space Flight Center ; Greenbelt , Maryland. 


located at 0.6 solar radii above the photosphere. This 
distance was chosen from a parameteric fit of this 
quantity based on comparisons of the model with the 
observed interplanetary field. This would be the location 
in the model where the highest coronal loops would 
form. Bugoslavaskaya [1950] observed the solar corona 
from 1887 to 1945, and Newkirk [1967] found the 
highest closed arches have a mean height of 0.6 solar 
radii above the limb. 

Further evidence for the highest closed magnetic 
loops lying near 0.6 solar radii above the limb is pro- 
vided by the observations of Takakura [1966] that U 
bursts have a maximum height near this value. U bursts 
are thought to be essentially type III radio bursts caused 
by the motion of high speed particles through the solar 
atmosphere, in which an increase in radio frequency 
emitted follows the usual decrease. The inversion in 
radio frequency emitted is interpreted as a decrease in 
altitude of these particles as they move through the 
corona on the magnetic field lines which govern their 
motion. 

Although the magnetic models of the corona of 
Altschuler and Newkirk [1969] and Schatten et al 
[1969] appear to be capable of calculating the large- 
scale structure of the coronal and interplanetary mag- 
netic fields moderately well, there are two areas 
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where notable deviations may be found that relate to the 
magnetic models: 



Figure 1 . Magnetic field geometry in the “source sur- 
face ” or “zero-potential surface ” models. The fields are 
constrained to the radial direction by the solar wind in 
these models . The equations obeyed in the different 
regions are shown . 

1, Solar flares appear to affect the large-scale magnetic 
field of the corona. The influence may appear in a 
solar eclipse photograph as the formation of series of 
fine rays directed radially away from the source of 
the flare [Smith and Schatten, 1970] . 

2. Although much of the open field structures and 
closed field structures have the correct topology, the 
structures are not always directed properly. A notable 
example is the polar plumes, which appear to bend 
continually equatorward, whereas the magnetic 
models orient them in the radial direction at the 
“source surface” or “zero potential surface”. Another 
example are streamers, whose axes show a prefer- 
ential lean to the equator near solar minimum and 
toward the poles at solar maximum [Waldmeier, 
1970]. Figure 2 illustrates the nonradial aspects or 
coronal features. 

The first area of disagreement is expected due to the 
large amount of hot plasma emitted by a flare. The 
current-free assumption in the inner corona is violated 
by this hot plasma and thus the potential solution is no 
longer valid. The second area of disagreement may relate 
to the latitudinal and azimuthal magnetic pressure terms 
that are important in coronal structure but have been 
neglected beyond the zero potential surface in prior 
work for mathematical simplicity. The purpose of this 
work will be to improve the model by including this 
effect. The energy density of the radial magnetic field 




Figure 2. Drawing of the February 15, 1961 solar 
eclipse (top) and the February 25, 1952 solar eclipse 
(bottom) superposed with a “source surface ” at 1.6 
solar radii. Npte that most closed arches fall below this 
surface . Beyond this distance most structures are “open” 
but not strictly radially oriented. 

(providing transverse pressure stresses) falls off much less 
rapidly than that of the transverse field [Schatten et al 
1969] . The energy density of the transverse field ap- 
proximately equals that of the plasma at about 0.6 solar 
radii (above the photosphere). Thus the plasma extends 
the magnetic field outward near this point. In the case of 
the radial field, quality with the plasma energy density is 
only reached at the Alfve'n point near 25 solar radii. 
Thus transverse magnetic pressure is expected to be an 
important effect long after the coronal plasma has be- 
come supersonic. The magnetic field behaves like open 
rigid wires along which the plasma is constrained to 
flow. The magnetic field thus may still guide the plasma 
motion from 0.6 to 25 solar radii. This paper suggests a 
method to mathematically calculate the magnetic struc- 
ture in this region. 
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CURRENT SHEET MODEL 

Out to the Alfven point, the value of p (ratio of plasma 
to field energy density) for the coronal plasma is signifi- 
cantly less than one. Thus currents flowing in the 
coronal plasma cannot apply a significant pressure on 
the magnetic field except where the field is weak (near 
regions of opposite polarity fields). Currents are neces- 
sary to open the magnetic field into sector-like struc- 
tures. If any significant transverse currents were located 
in regions of moderate field strength, a strong jXB force 
would occur which the plasma could not resist. Thus, 
the currents tend to be present in high p regions, where 
the field reverses, and the jXB forces are small. For these 
reasons, the transverse currents flowing in the coronal 
plasma in this model are constrained to flow only where 
the field is weak (near zero), hence on sheets near 
oppositely directed field regions. 

The magnetic model can be used for computing the 
location and strength of these current sheets as follows. 
The magnetic field is first calculated directly from 
potential theory by means of Legendre polynomial tech- 
niques to a particular surface — for example, a sphere of 
radius 1.6 solar radii. Although not utilized in exactly 
the same way, this sphere will again be referred to as the 
“source surface.” The magnetic field then is reoriented 



Figure 3. First step in the current sheet magnetic 
model A potential solution is derived for the field be- 
tween the “ source surface” and the photosphere. The 
field computed on the source surface is then reoriented 
so that it points outward everywhere. The field is then 
computed beyond the “source surface” from potential 
theory. The sense of the magnetic field is opposite half 
the time to what it should be. This “error” is corrected 
in the next step. 


such that it points outward everywhere; however, it is 
still along the same direction and possesses the same 
field magnitude. Thus if B r > 0 on this “source surface” 
the field is unchanged, but if B r < 0 then B r , Bq , and 
B<$ are replaced by -B r ,-BQ, and -Bq, The field is then 
calculated beyond the “source surface” from potential 
theory, again using a Legendre polynomial expansion of 
the field (see appendix). Now the monopole term is non- 
zero and rather large; thus it appears as if the sun has a 
high magnetic monopole moment and all the Legendre 
polynomial coefficients bear little or no relationship to 
their previous values (fig. 3). The physical effect is that 
beyond the “source surface” the magnetic fields cannot 
now form closed arches as they are all directed outward. 
This temporarily violates V • B = 0 on the source surface 
but this error will be corrected in a later step. This 
change of field direction does not affect the magnetic 
stresses: They will remain the same across the “source 
surface” and the field will still form a minimum energy 
configuration (with the condition that the field lines 
remain open). The last step (fig. 4) is to return the 



Figure 4. Second step in the current sheet magnetic 
model The field that was disoriented is reoriented by 
reversing the sense of the magnetic field components. 
This requires a current sheet to be employed in the 
corona to separate regions of oppositely directed field to 
obey Maxwell's equations. Allow the magnetic field to 
“open” by thin current sheets is consistent with the 
physical model of this region of the corona possessing a 
low p. If significant transverse currents flowed elsewhere 
a strong jXB force would develop which the plasma 
could not maintain. This model may be used to calculate 
the magnetic oriented structures in the corona with less 
simplified solar wind currents . 
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magnetic field to its former sense of direction with the 
calculated strength and orientation. This violates 
VX B = 47rj jC- 0 unless appropriate current sheets are 
introduced as shown. Physically, current sheets are intro- 
duced between areas of oppositely directed fields and 
thus prevent the field from forming arches beyond the 
“source surface.” Note that the polar fields and streamer 
fields possess similar shapes to those in the corona 
(fig. 2) and not the radial orientation seen in figure 1 . 

The invariance of the Maxwell stress tensor under this 
field reversal scheme is important to ensure against 
unequal stresses across the “source surface.” The 
Maxwell stress tensor is defined such that j X B = V • M , 
The stress tensor is shown in figure 5 . As can be seen, 
changing the sign of the three components leaves M 
unchanged; thus, the magnetic stresses in the corona are 
balanced. 

T x B = V -to 
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Figure 5. The Maxwell stress tensor . Note that it is 
identical if all three components are reversed. This 
allows the stresses to be balanced after the field reversal 
processes. 


COMPARISONS OF THE CURRENT SHEET MODEL 
WITH OTHER MODELS AND THE SOLAR CORONA 

The current sheet model is first compared with the 
“source surface” and the “zero potential surface” 
models as well as an exact MHD solution for an axisym- 
metric isothermal corona. This latter solution has been 
computed after the formalism of Pneuman and Kopp 
[1970] for the corona with a temperature of 
1.56X10 6 °K and a dipole field. Figure 6 shows this 
comparison with an assumed dipolar solar field. The 
field lines labeled with crosses represent the present 
study with the “source surface” located at 1.6 solar 
radii. Solid lines indicate field directed away from the 
sun, and dashed lines indicate field toward the sun. The 
heavy solid lines indicate the MHD isothermal coronal 
solution of Pneuman and Kopp. The dashed and dotted 
lines indicate the field lines calculated by the Altschuler 
and Newkirk model with a zero potential surface located 
at 2.5 solar radii. The “source surface” solution of 
Schatten etal. [1969] is similar to this solution except 
the field lines would be oriented radially somewhat 
closer to the sun. As can be seen, the field lines com- 
puted from the isothermal MHD solution and the cur- 
rent sheet solution are nearly identical. The foot points 


Figure 6. A comparison of solutions for a dipole solar 
field with the zero potential surface model [Altschuler 
and Newkirk , 1969] the present model and the exact 
isothermal MHD coronal solution [Pneuman and Kopp , 
1970]. Note that the present solution is quite similar to 
the isothermal coronal solution supporting the suggestion 
that the coronal currents are confined to thin sheets. 

of the field lines indicate the quality of their agreement. 
The magnetic potential solution begins to diverge from 
the other solutions near the zero potential surface. The 
rather close agreement between the current sheet solu- 
tion and the MHD solution suggests that much of the 
current flowing in an isothermal corona does so near 
current sheets as suggested earlier. Altschuler and 
Newkirk [1969] chose the location of the zero potential 
sphere to be 2.5 solar radii based on a comparison of 
field geometry with coronal forms, whereas Schatten et 
al. [1969] chose the 1.6 solar radii value for the “source 
surface” based upon the observed highest closed arches 
and agreement with comparisons of their model with the 
interplanetary magnetic field. In the present model, if 
the “source surface” is set at 1.6 solar radii, the shapes 
of features are similar (out to 2+ solar radii) to the 
Altschuler and Newkirk result, and the coronal magnetic 
field extends out from 1.6 solar radii, similar to the 
result of Schatten et al [1969] . Thus the disagreement 
between these two values where the coronal magnetic 
fields extends outward may be ended by utilizing this 
new model. The agreement with the axisymmetric MHD 
solution suggests that the current sheet model may now 
be used with more confidence in calculating fields in 
three dimensional nonsymmetric situations as well. 

First, however, let us examine whether the current 
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sheet model can calculate nonradial streamers and com- polar plumes. The field pattern beneath shows calcula- 

pare them with observed nonradial streamers. Figure 7 is tions from the current sheet model using an axisym- 

a drawing from Bohlin of the May 30, 1965, solar eclipse metric magnetic condition. As can be seen, rather non- 

from photographs by Smith (top). This eclipse shows radial field lines may be computed in the model quite 

several nonradial streamers in addition to the nonradial similar in appearance to the structures observed. The 
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Figure 7. A comparison of the structure of the solar corona during the May 30, 1965 
solar eclipse (top) with computations from an assumed axisymmetric photospheric 
field pattern (bottom). The shape of the polar plumes is calculated quite well in this 
model as well as a nonradial helmet streamer. 
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polar field lines appear similar to the polar plumes. The 
computed field configuration in the equatorial regions 
are also oriented toward the equator as in the eclipse 
drawing. 

A computation of the magnetic field projected into 
the plane of the sky from this model for the November 
12 , 1966, solar eclipse is shown in figure 8 superimposed 
with a drawing of the coronal forms by Newkirk et al 
[1970] , The solid lines indicate away-from-the-sun mag- 
netic field, and dashed lines indicate toward-the-sun 
field. Many of the features line up surprisingly well with 


the field lines calculated as can be seen. Large-scale mag- 
netic loops are calculated near streamers j3 ' and f and 
closed arches are observed underlying these streamers. 
Many of the “open” magnetic field lines near regions j3, 
f, and 7 are closely alined with coronal features in the 
same areas. The general agreement of the magnetic field 
calculations with the observed features for this solar 
eclipse is rather good. 

CONCLUSIONS 

A new current sheet magnetic model for the solar corona 



Figure 8. A comparison of the computed coronal magnetic field with a drawing of the solar eclipse features from 
Newkirk and Altschuler [1970] for the November 12, 1966 solar eclipse . The solid lines represent magnetic field 
directed away from the sun and the dashed lines field toward the sun . The field lines originate on the photosphere 
in the center of each of the 648 grid points. Many of the observed features line up with the computed field lines. 
The computed field lines terminate at five solar radii and are projected into the plane of the sky. 






has been developed. It is capable of calculating the quiet 
large-scale magnetic field structure in the corona. As sug- 
gested by physical arguments, thin current sheets are 
utilized to separate regions of oppositely directed fields. 
This approximation appears to be a rather good one in 
that the dipole solution is nearly identical with the iso- 
thermal MHD coronal solution of Pneuman and Kopp 
[1970]. 


A comparison of field computations with the ob- 
served structure for the May 30, 1965, solar eclipse re- 
veals that the model appears to be capable of calculating 
the orientation of the polar plumes fairly well, as well as 
nonradial streamer configurations. A comparison be- 
tween the computed magnetic field and the observed 
solar corona for the November 12, 1966, solar eclipse is 
shown. Many of the observed features are also seen in 
the computed magnetic field. 


APPENDIX This appendix discusses the solution of the field beyond the “source surface” method of 
fitting the vector Legendre polynomial coefficients to the three dimensional vector field 
on the “source surface.” The vector field up to and including the “source surface” is 
computed in accord with the techniques of Altschuler and Newkirk [1969] without using 
any currents in the solar corona. The present model may be improved in the future by an 
iterative process using the currents computed in the present model to calculate the 
solution below the “source surface” as well and then recomputing the field beyond the 
“source surface.” This may represent a minimum improvement, however. 

The solution for the Laplacian equation in spherical coordinates for r > R is 

A 

Wr,8,4>) = R ZX {(T) [*«* cos(mtf>) + ^ s in(m0)]/^(0)J (1) 

n=o m=o 


The components of the magnetic field are: 



n - L M 

n< t > r sin 0 d<p 


-/•fcy.o 


(4) 


The associated Legendre polynomials utilizing the Schmidt normalization have been 
used [Chapman and Bartels , 1940] . Thus to determine the magnetic field beyond the 
“source surface” it is necessary to compute ^ and h™ from the vector field on the 
“source surface” as a boundary condition. 
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The components of the magnetic field on the “source surface” are first oriented 
away-from-the sun so that if B r < 0 on the “source surface,” the signs of B r , Bq, andB^ 
are reversed. 

In this analysis we have utilized a photospheric grid of 27 longitudes and 24 latitudes 
in equal steps of sine (latitude). We have also chosen N— 9 as the maximum principal 
Legendre index to consider for practical considerations. A least-mean-square fit to an 
overdetermined linear system of 1944 (27X24X3) equations involving 100 unknowns is 
then utilized to best fit the vector field on the “source surface,” we let 



where B(i,j,k) equal the vector field components, and k = 1, 2, or 3 refers to the radial, 
latitudinal, or azimuthal field component at Of and <f>j. 

It is necessary now to obtain the ^ and h™ that minimize F\ the sums of squares of 
the differences between the known components of the field on the “source surface” 
B r (iJ), BqQJ) and B^QJ); and the component values computed from ^ and h ™ at Q[ 
and <pj. 

Let us choose 


a ntm = ( w + 7 ) cos m <t> P%( e ) 

$nm\~ (n + 7 ) sin m( t> P%( e ) 

a nm2 = - c ™ m< l> £ P%( 9 ) 

> 

a nm3 = -£r e Sinm<t> ^ (0) 
cos m<l>F%(e) 

Thus equation (5) becomes 

-EEE B(i, j, k ) - ^7 cmmk + h™ finmlc) 

i j k L n m 


( 6 ) 


( 7 ) 


The equations to minimize F are: 


-^- = 0^ = 0 
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For each («, m) may be rewritten: 


/ / l r=o s-o L J J (9) 

EEE 1^(4 /, ^)P nm k^y^nmk^ 2 + ® 

i j h \ t=o s=o ' f (10) 

where fands are dummy indices used for n andm. The unknowns areg^ 2 and and 
B(iJ,k) is the known vector field; a and J3 are known from equation (6). 

The column vectors 


%<p2 1 


001 01 1 
00102 1 


“001 

“ 010,2 


= %01 

^ 0 , 0,2 


are of length 24X27X3 for each 7VM and 


5 ( 01 , 01 , 7 ) 


5(0 ,,0 2 ,/) 


5(0, 0,7) 
5(0i,0i,2) 


5(0,0, 5) 
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Now defining the matrix a0 such that the rows of o0 are as follows: 


oj3(l)= a 00 
00 ( 2 ) = a i0 
a0(3) = an 


o^(55) = o g9 
a0(56) -0i i 

00 ( 100 ) 

with all m = 0 elements missing from h™ and from 0; GH is a 100X 1 matrix and a0 is a 
100X 1944 matrix. Equations (9) and (10) may be rewritten as 


A0 ■ B =AB ' GH (12) 

choosing AB(iJ) ~ cx0(i) • o0(f) so that o0 is a 100X 1944 matrix, B is a 1944X 1 matrix, 
ABfj is a 100X 100 matrix and GH is a 100X 1 matrix. 

By an inversion of the symmetric matrix AB, GH may be solved as follows: 

GH - AB _1 * o0 * B 

This requires inverting a square matrix each of whose sides equals (N + l) 2 which for 
N = 9 is 100, yielding estimates for $ and h™ that arise from a least mean square fit to 
the three vector components of the magnetic field on the source surface. Equations (2), 
(3), and (4) allow a computation of the magnetic field everywhere above the “source 
surface.” It is necessary, however, to reverse the sense of the three components of the 
magnetic field depending on whether the footpoint of the field line has had its sense 
reversed (if B r < 0). Those field lines are shown as dashed lines in the figures. 
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DISCUSSION R. A. Kopp I just wanted to know how you got the field lines to go through the 
source surface radially. You said you did that calculation first with the ordinary source 
surface. Don’t the field lines have to go perpendicular as they pass through the source 
surface? 

K H. Schatten No, they weren’t taken to be radial in this case. It was strictly a 
potential theory using the Legendre polynomial coefficients, and we did not constrain 
them to be radial - they just occur in whatever direction they happen to be there. Then 
we did a least-mean-square solution so that the three vectors, B r B® B^ ; are at the source 
surface, so that there was no necessary restriction into the radial direction there. 

R. B. Leighton It’s a kind of heuristic way of recognizing you’re going from a case 
Where there is essentially a potential magnetic field problem to a fluid problem but in 
which the magnetic field still plays a role in terms of its stress, in guiding the flow. This is 
why you get away with reversing the field there and keeping on going and damn the 
torpedoes — right? 

K. K Schatten Right. 

F. C. Michel I don’t know if it’s directed to Ken Schatten or to Gordon Newkirk. But 
the eclipse of 1970 looked quantatively as if it had been drawn by a student who hadn’t 
learned his lessons very well, and if he had drawn it in the form it appeared he would be 
critized for having a rather poor representation of how the lines should go. There are great 
regions where the striations are essentially parallel; and not only that but they didn’t 
appear at all radially from the sun, and there were about five such regions. I just won- 
dered, what’s your attitude towards this? Do you think this is a significant deviation from 
the theory or do you think that’s easy to account for? 

K . H. Schatten Yes, I think it is a significant deviation. As Gordon Newkirk men- 
tioned, there was a lot of solar activity at the time. These models essentially assume 
quasistationary conditions with no sort of solar activity to mess things up. When you have 
activity you get a lot of plasma ejected into the corona, and currents can form which will 
twist the field to almost any configuration, and you can’t just compute it from the 
photospheric field. You can imagine a very strong flare occurring in a particular location 
where you have, say, a stream or something like that, and it blows all the field lines out 
into radially away from that point, or into some other configuration. That probably is 
why the solar eclipse calculations for 1970 were not as good as some of the others. 
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